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It is a half-century-long story about the first-order phase transition effects 
on equilibrium, stability and pulsations of planets and stars. The topics more 
or less considered in author's papers are mainly touched, and no attempt was 
q\ • done to cover all literature on subject. 

ON , It was Ramsey who in 1950 had first shown in MNRAS paper [1] that if in 

the center of a planet, with central pressure P c just reaching critical pressure 
Feb the first-order phase transition takes place with density jump from pi to 
P2 = Q* Pi, and if q > 1.5, then planet loses its stability at the same moment. 
This remarkably amazing and quite unexpected result has been since then 



rediscovered by many authors repeatedly - last time in nucl-th /9902033 
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j> ■ D'you think all these (cloths) will be the wear? 

ON , - / think all these should be tailored. " 

q \ Yuri Levitanski, Soviet Jewish poet 



The story began to me some 35 years ago when trying to understand paper [2], I've 
found that something is wrong with Mass-Central Density curves for white dwarf stars 
near maximum - there should be sharp, not smooth maximum of mass, as reverse /3-decay 
Q-f reactions lead to the discontinuity of density distribution inside the star. 

To simplify the problem, I used some models, allowing analytical investigation - 
polytropes with indices n equal to and 1 in envelope and core in various combinations. In 
all cases considered, it happened that if q > 1.5, the instability occurs at the moment the 
phase transition starts in the center of the star. 

My supervisor acad. Ya.B. Zeldovich for a long time did not believe in such "strange" 
figure 3/2, so I kept analyzing various models. After a while Ya.B. became himself convinced 
and within a brief period he developed an amazingly elegant method of prooving this word 
constant 3/2, and a joint paper [4] was submitted to Astronom. Zhurnal. Needless to say, 
I was happy that I was right and that I had joint paper with such an outstanding scientist. 
Catastrophe came short: Ya.B. sent me a brief letter noting that "the result 3/2 is known 
in literature!". 
No comments... 

I've never managed since then to write another joint paper with Ya.B.... 
Why 3/2? I do not know, but it seems that 3/2 number is directly related to the r _1 - law 
of potential U(r) in Newtonian Theory of Gravitation (NTG) in 3D space[ll]. 
This is not the last surprise of Phase- Transition Astrophysics (PTA)! 
In the classical theory [5] of equilibrium and stability of stars comprised of matter with 
smooth Equations of State (EoS), there is a common sense that rotation leads to the increase 
of stability reserve while General Relativity (GR) leads to the decrease of stability. 

And you may guess that exactly opposite situation is in PTA. In GR, the critical value 
of energy density q = e 2 /£i is [6]: 3/2 (1 + P /ei), that is larger than in NTG. 
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Why larger? - I do not know. . . 

As to rotation, it was found [7], that for steady-state rotation with small angular velocity 
Qcrit = 3/2 — Q 2 I 4irGpi, that is, of course, rotation reduces the stability of star against 
phase-transition- induced instability. 
Why reduces? - I do not know... 

Abovementioned three suprises of PTA, combined in the formula 

q crit = 3/2 - n 2 I AnG Pl +3/2(1+ P / P iC 2 ), 

are valid for any EoS's of old and new phases. 

A number of amazing results was found analysing various particular models. First to 
be mentioned is the two-constant-density-phase model with First-Order Phase Transition 
(PT1) at the boundary between core and envelope. This model was extensively used for 
different problems - general dependence of Mass-Radius etc. curves, effects of rotation and 
GR, neutral core (when PT1 begins at some distance from the center of a star) etc. 

In the last case, it was found [3] that q crit is an increasing function of size of neutral 
core: q crit = 1 + fc/[3 — 4 x — (k — 1) a; 4 ], where x is relative radius of neutral core, and k 
is relation of density in neutral core to density in envelope; and q crit — > oo at x — > x cr (k), 
where for example at k — 1, x cr = 3/4 - another amazing value. At larger neutral cores, 
PT1 with arbitrary large q can not force a star to lose its dignity and stability! 

For a model with polytropic indices n — 1 both in envelope and neutral core [3], x cr = 
.6824 (corresponding relative mass of core is .6375). 

Returning to n = 0, for q > 3/2 there is another critical point in Mass — P c curves, 
where at the minimum of mass the recovery of stability takes place, and for larger P c there 
is a branch of stable equilibrium states. At theminimum of Mass: 

f(q,x) = (q-l) 2 x 4 + 4(g-l)x + (3-2g) = 0. (1) 

It was found [8] that GR effects, in the first Post — Newtonian approximation lead to the 
PN-correction to x in Eq. (1): 

9 - 7 q + 27(g - l)x + (q - l)(4g - 27):r 2 + (g - 1)(9 - 4q)x 3 P 
PN{q ' X) ~ 2(g-l)[l + (g-l)*3]3 — 2 ■ ( 2 ) 

The surprises have not finished - this correction is negative at small values of q and positive 
at q > 1.89. That is for larger q, GR effect is of correct sign and reduces the region of 
stability at (q — x) plane. 

For the same n = model, the analytical formula for frequency, u, of small adiabatic 
radial pulsations of the lowest mode can be found: 

Wo " 3( g -l)(l-s)' (3) 
and in the case of slow rotation with angular velocity Q [10]: 

= w o + A n (q,x), with 



An(?,a;) = |n 2 



5x(l - x)(l + x) 2 l + {q-l)x 



l + (q-l)x 5 (q-l)(l-x) 



(4) 
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Rotational correction to frequency squared is negative at smaller values of x, that is 
rotation reduces the stability reserve of star with PT1, then for larger value of x, rotational 
effect is of "correct" sign. 

In fact, dependence of both corrections, due to GR and rotation, on x and q is rather 
complicated, and Fig. 1 presents only a part of (q—x) — plane with lines on which A^(q, x) — 
(broken line labeled as "Rot")and A PN (q,x) = (dash-point line labeled as "GR"). Also 
shown is the curve f(q,x) = (solid line labeled as "Crit") from Eq.l, which marks a 
boundary between stable equilibrium states (right-hand) from unstable ones (left-hand). 
Remarkably, this curve crosses both curves of zero correction. Due to GR, curve f(q, x) — 
is forced to rotate clockwise around point of intersection of lines " GR" and " Crit" , while 
rotation makes the critical curve rotate counterclockwise around point of intersection of lines 
"Rot" and "Crit". 
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FIG. 1. Three important curves at (q — x)— plane: critical state of star stability recovery, see 
Eq. 1, and GR and Rotational corrections equal to zero, see Eq. 2 and Eq. 4, respectively. 

Why do rotation and GR exert such strange effects on stability of a star with PT1? 

- I do not know. . . 

A lot of interesting and unsolved problems of pulsations (eigen- values and functions, 
damping), strong rotation with regard to deviation of equilibrium figures from spherical 
symmetry, considering more EoS, etc. have been left aside , but 1500-word mark is close 
and I pass to the epilogue... 

Most of this story happened to me many years ago, in the Soviet Union, the Power so 
unexpectedly crushed in the recent past, as if some first-order phase transition had acted in 
the center of It... 

Now I'm in Israel, all my papers being left over there and these lines being written by 
heart, by memory, with no paper at hand... 

And last sentences (hopefully still in 1500- word limit), returning to epigraph: 

- Do I believe this all will be awarded? 

- I think this all should be written. 
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